arXiv:astro-ph/0211263vl 12 Nov 2002 


Astronomy Astrophysics manuscript no. borkova February 2, 2008 

(DOI: will be inserted by hand later) 


Two populations among the metal-poor field RR Lyrae stars 

Borkova T.V., Marsakov V.A. 

Space Research Department at Rostov State University, 

and Issac Newton Institute of Chili Rostov-on-Don Branch Stachki 194, Rostov-on-Don, 344090, Russia 
e-mail: borkova@rsusul.rnd.runnet.ru, marsakov@ip.rsu.ru 

received 17 April 2002 ; accepted 17 September 2002 

Abstract. We compute the spatial velocity components and the galactic orbital elements for 209 metal-poor 
([Fe/H] < —1.0) RRLyrae (ab) variable stars in the solar neighborhood using proper motions, radial velocities, 
and photometric distances available in the literature. The computed orbital elements and published heavy element 
abundances are used to study relationships between the chemical, spatial, and kinematical characteristics of nearby 
field RR Lyrae variables. We observe abrupt changes in the stellar spatial and kinematical characteristics when the 
peculiar velocities relative to the local standard of rest cross the threshold value, Upec ~ 280 kms^ 1 . This provides 
evidence that the general population of metal-poor RRLyrae stars is not uniform, and includes two spherical 
subsystems occupying different volumes in the Galaxy. Based on the agreement between typical parameters of 
corresponding subsystems of field RRLyrae stars and of the globular clusters, studied by us earlier, we conclude 
that metal-poor stars and globular clusters can be subdivided into two populations, but using different criteria 
for stars and clusters. We suppose that field stars with velocities below the threshold value and clusters with 
extremely blue horizontal branches form the spherical, slowly rotating subsystem of the proto-disk halo (related 
by its origin to the Galactic thick disk). It has a negligible, but non-zero, vertical metallicity gradient. Field stars 
with fast motion and clusters with redder horizontal branches constitute the spheroidal subsystem of the accreted 
outer halo, with is approximately two times larger in size than the first subsystem. It has absolutely no metallicity 
gradients, most of its stars have eccentric orbits, many stars display retrograde motion in the Galaxy, and their 
ages are comparatively low, supporting the hypothesis that the objects in this subsystem have an extragalactic 
origin. 
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1. Introduction 

The presence of two different populations with separate 
histories in the metal-poor halo was suggested by Hartwick 
(1987). He showed that model of the dynamics of RR 
Lyrae variables with metallicities [Fe/H\ < —1.0 required 
two components: one spherical and a somewhat flattened 
component that is dominant at galactocentric distances 
less than the radius of the solar circle. The idea that 
there are two subsystems in the metal-poor halo has al¬ 
ready been addressed in several investigations focused on 
globular clusters. It turned out, that globular clusters 
present a distinctive inner feature (the morphology of their 
horizontal branches), which makes it possible to distin¬ 
guish the individual metal-poor clusters in terms of dif¬ 
ferent halo subsystems. It was found that halo clusters, 
that have redder horizontal branches for a given metal¬ 
licity (i.e. with horizontal branches showing a consider¬ 
able number of stars on the low-temperature side of the 
Schwarzschild gap), are predominantly outside the solar 
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circle. Furthermore, they exhibit larger velocity disper¬ 
sion, slower circular velocity (a significant number hav¬ 
ing retrograde orbits), and are, on average, younger than 
clusters with extremely blue horizontal branches, which 
are concentrated within the solar circle (Da Costa & 
Armandroff 1995; Borkova & Marsakov 2000). Note that 
relative ages for globular clusters were determined by dif¬ 
ferent authors with the help of observed high-precision 
colour-magnitude diagrams and theoretical isochrones. 
Therefore, these ages are independent of horisontal branch 
morphology, and the relative youth of some metal-poor 
globular clusters is beyond doubt (see compilative cata¬ 
logue of homogeneous age dating of 63 globular clusters 
by Borkova & Marsakov (2000)). The explanation sug¬ 
gested for the difference between these two populations 
was that the subsystem of clusters with extremely blue 
horizontal branches (i.e. the older halo) formed together 
with the Galaxy as a whole, whereas the clusters of the 
younger halo subsystem formed from fragments captured 
by the Galaxy from intergalactic space at later stages of 
its evolution (Zinn 1993). Recent observations (Ivezic et 
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al. 2000, Vivas et.al. 2001) suggest strong evidence in fa¬ 
vor of the hypothesis that star and also globular clusters 
of the outer halo are the debris left over from the accretion 
of dwarf galaxies. Unfortunately, it is impossible now to 
ascribe concrete nearest field RR Lyrae stars to some sub¬ 
system, because we do not know any intrinsic distinctive 
quality for them. 

Some papers (Chiba & Yoshii 1998; Martin & Morrison 
1998; Dambis & Rastorguev 2001) present detailed studies 
of the kinematics of RR Lyrae stars in the solar neighbor¬ 
hood. The stars were assumed to form only one subsys¬ 
tem in a metal-poor halo. We use spatial velocities and 
computed elements of galactic orbits as our main criteria 
for isolating subsystems (because of the local, near-solar 
position of the studied RR Lyrae stars). Our wish is to 
investigate relationships between physical, chemical, spa¬ 
tial and kinematical characteristics of RR Lyraes stars in 
each metal-poor subsystem, determine the characteristic 
parameters of these subsystems, and compare them with 
the parameters of similar subsystems of globular clusters. 

2. The Data 

We used for this study the largest catalogue of RR Lyrae 
variables, compiled by Dambis & Rastorguev (2001). The 
catalogue contains 262 stars with published photoelec¬ 
tric photometry, metallicities, radial velocities, and abso¬ 
lute proper motions. Dambis & Rastorguev (2001) used 
the proper motions from Hipparcos, PPM, NPMI, and 
Four-Million Star Catalog (Volchkov et al. 1992); the 
proper motions from the last three ground-based cata¬ 
logs have been reduced to the Hipparcos system. The 
metallicities were used mainly from Layden (1994) and 
Layden et al. (1996), whose metallicity scale is in con¬ 
formity with metallicity scale of globular clusters (Zinn 
& West 1984). Dambis and Rastorguev, (2001) used ra¬ 
dial velocities mainly from the paper of Fernley et al. 
(1998), and Solano et al. (1997), and main magnitude 
< V > - mainly from Fernley et.al. (1998). Bearing 
in mind the large relative errors of trigonometric paral¬ 
laxes for distant objects, we chose the photometric dis¬ 
tance scale of Dambis & Rastorguev (2001), assuming 
M V (RR) = 0.26 [Fe/H] + 1.17. For each star, we com¬ 
puted the spatial velocity components in the cylindrical 
coordinates and the orbital elements using the Galaxy 
model from Allen & Santillan (1991), which includes a 
spherical bulge, a disk, and an extended massive halo. The 
model assumes the galactocentric distance of the Sun to 
be Rq = 8.5 kpc and the Galactic circular rotation ve¬ 
locity at the solar distance to be V I0 t = 220 kms -1 . The 
orbital elements were computed by modeling five com¬ 
plete revolutions around the galactic center for each star. 
The most informative quantities are: Z max - the maxi¬ 
mum height of the star above the galactic plane, R a - the 
orbital apogalactic radius; R p - the orbital perigalactic ra¬ 
dius, and the eccentricity, e = ( R a — R p )/(R a + R p ). A 
choice of stars based solely on their variability type and 
their visible magnitude ensures an absence of kinematical 


selection effects in the catalogue used. Our final sample 
of the metal-poor RR Lyrae variables contains 209 stars 
with metallicity [Fe/H] < —1.0. This simple metallicity 
criterion eliminates objects that belong to the thick disk 
subsystem of the Galaxy (see substantiations of Borkova 
& Marsakov (2002)). 

3. Criteria to select the subsystems 

It is much more difficult to identify objects that have an 
extragalactic origin, i.e. those belonging to the accreted 
halo. According to the hypothesis that the protogalaxy 
collapsed monotonically from the halo to the disk, sug¬ 
gested by Eggen et al. (1962), stars genetically related to 
the Galaxy cannot have retrograde orbits. Only the old¬ 
est halo stars may be an exception, since they could have 
retrograde orbits due to the natural initial velocity dis¬ 
persion of protostellar clouds. On the other hand, some 
stars formed from extragalactic fragments and captured 
by the Galaxy may have a prograde orbits. In any case, 
such stars should have fairly large peculiar spatial veloc¬ 
ity relative to the local standard of rest, U pec - Fig. |ja dis¬ 
plays the relation between the peculiar velocity (assuming 
( U 0 , V 0 , IF o )=(-10,10,6) kms -1 and the azimuth (circular) 
velocity component, O, for the RR Lyrae stars of our sam¬ 
ple. The diagram shows that there is a transition from pro¬ 
grade to retrograde orbits around the galactic center near 
Upec ~ 280 kms -1 . We also observe an abrupt increase 
of the dispersion of the circular velocity component and a 
break of the dependence of circular velocity as a function 
of peculiar velocity at the same place (see the dispersion 
bars and the regression lines in the diagram). Fig. [ijb dis¬ 
plays a significant increase in the scatter of stars in Z max 
when crossing the same threshold peculiar velocity. The 
abrupt change in the apogalactic radii of the stellar or¬ 
bits is even more evident (Fig. [l]c). First, the mean value 
of Ra remains much the same and when peculiar velocity 
crosses the threshold level, the mean value of apogalactic 
distances monotonically increase and their scatter sharply 
extends. We can see in Fig. [Tel that the perigalactic dis¬ 
tances also demonstrate the break of their relations in the 
vicinity of the same threshold peculiar velocity. The or¬ 
bital eccentricities (Fig. ||c) not only abruptly increase 
their dispersion when crossing the same point, they also 
demonstrate different relations. First, the orbital eccen¬ 
tricities increase almost linearly with the peculiar veloci¬ 
ties, reaching a maximum near the threshold velocity level. 
With further increase of Vj, e c, the mean and the scatter 
of the eccentricities do not change, within errors. Fig. [If 
shows that RR Lyrae stars with velocities near the thresh¬ 
old level might have any orbital inclination up to orthogo¬ 
nal to the galactic plane, whereas the range of ” permitted” 
inclinations continuously becomes narrower when velocity 
moves away from this threshold level on both sides. For 
these reasons, we adopt Vp ec > 280 kms -1 as the critical 
value for distinguishing stars of the outer accreted halo. 
Here we suppose that the stars with lower peculiar ve¬ 
locities have a galactic origin, and belong to the proto- 
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Fig. 1. The relations between the peculiar velocity relative to the local standard of rest, V pec , and other characteristics 
of the RR Lyrae variables. Large open circles with error bars are mean values and dispersions of the corresponding 
parameters in narrow intervals y pec . The straight lines are the least-square fits for the stars lying both to the left and 
on the right of V pec « 280 kms -1 (the vertical dotted lines). An abrupt change in the relations and the dispersions 
near this threshold level is apparent in all panels. 
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disk halo subsystem. Apparently, this kinematical crite¬ 
rion is not entirely unabbiguous: some stars of the proto¬ 
disk halo may have larger residual velocities. Evidence for 
this is provided, in particular, by the increase in the stellar 
density immediately to the right of the thereshold veloc¬ 
ity level in our diagrams. However, we decided to retain 
a simple criterion, in order not to artificially confuse the 
situation. 

Note again that the principal criterion distinguishing 
globular clusters of the accreted halo is their redder hor¬ 
izontal branches compared to clusters of the proto-disk 
halo. Moreover, empirical evidence suggests that the mor¬ 
phological structure of the horizontal branches is related 
to the mean period of RR Lyrae variables in metal-poor 
clusters (vanAlbada & Baker 1972). We investigate the de¬ 
tailed data for RR Lyrae in globular clusters from Clement 
et al. (2001), and show that the redder the horizontal 
branch, the lower the mean period with a height corre¬ 
lation coefficient r = 0.8±0.1. It means that the mean 
period for RR Lyrae in clusters can be used to divide 
this cluster into subsystems. Therefore, it is interesting 
to test whether the variability period of single field RR 
Lyrae star may be adopted as well. To test this for field 
stars, we computed the mean periods of stars in a narrow 
metallicity range (—1-7 < [Fe/Hj < —1-2) for both halo 
subsystems. Each sample contained about sixty stars. It 
turns out that their mean periods are the same. Thus, the 
variability period of the field RR Lyrae stars cannot serve 
as an additional (internal) criterion to divide them into 
different subsystems of the nretal-poor halo. 

4. Properties of stars in the subsystems 

Let us now compare the properties of the resulting sub¬ 
systems. The spatial velocities of the stars can be used 
to obtain an estimate for a number of characteristics of 
subsystems, if we can first reconstruct the stars’ galactic 
orbits. Fig. |] shows the distributions of RR Lyrae in the 
two metal-poor subsystems as a function of their orbital 
elements. The top panels present the histograms of the ro¬ 
tation velocities for the sample of stars in the proto-disk 
halo and in the accreted halo. It is worth nothing that 
the two distributions can be properly fitted by a Gaussian 
(see the solid lines). Their maxima are separated by al¬ 
most the dispersion (Tabic |l|) . The second row of graphs 
in Fig. H shows the corresponding distributions in orbital 
eccentricities. Here the characters of the histograms are 
obviously very different. All eccentricities are present in 
approximately equal number in the proto-disk halo (only 
a small excess is seen towards the high eccentricity side). 
Stars with very eccentric orbits prevail in the outer halo 
(where almost two thirds of all stars have (e > 0.8)). The 
next row (Fig. [|c,f) presents distributions of orbital incli¬ 
nations. The stars of both subsystems can have any orbital 
inclination. In both cases, the number of stars strongly in¬ 
creases with decreasing inclination; however, this is true 
for stars with prograde orbits in the proto-disk halo, and 
for those with retrograde orbits in the outer halo. We must 


bear in mind that the deficiency of stars with large orbital 
inclinations is largely due to the kinematical selection ef¬ 
fect imposed on the sample of nearby stars. The vertical 
components of the spatial velocities of such stars in the 
solar neighborhood should be comparable to the galactic 
rotational velocity at this distance. Thus, the probability 
of their presence here is very low. 

The fourth row of the histograms (Fig. ||g,h) can be 
used to estimate radial sizes of the subsystems and the 
fifth row (Fig. ||i,j), their vertical sizes. To quantitatively 
estimate the outer sizes of the subsystems based on these 
distributions, we reject the five most distant points in 
each histogram. Therefore, we simultaneously remove the 
largest uncertainties in the determination of the orbital 
elements and avoid possible errors in assigning some stars 
in our sample to a particular subsystem. Such estimates 
indicate that the size of the outer halo is approximately a 
factor of two larger than that one of the proto-disk halo 
(Table [j]). The sizes of the subsystems in the direction 
perpendicular to the galactic plane also differ drastically 

(Table 0 ). 

It is obviously not correct to compute the scale height 
using Z max , since all stars of a subsystem cannot simulta¬ 
neously be located at the the highest points of their orbits. 
To reconstruct the ” real”, instantaneous Z-distribution for 
all stars, we must ’’spread” each star over its orbit from 
—Z m ax to +Z max in proportional to the probability den¬ 
sity of its location at different Z. (This operation also 
smoothes fluctuations in the histogram due to the limited 
number of stars in the sample.) This probability density 
can easily be found from the computed orbit of the stars. 
The details of the procedure were written by Marsakov & 
Shevelev (1995). The filled dots in the histograms of the 
bottom row (Fig. Bo) are the final reconstructed distri¬ 
butions in Z. In other words, this is how the stars will 
be distributed in height after some time, if they are ran¬ 
domly distributed in their orbits. The solid curves in Fig. 
|ji,j are least-square approximations of the reconstructed 
distributions using an exponential law, 

n = a ■ e _z / z °, 

where Zq is the scale height (the corresponding Zq val¬ 
ues are indicated in Fig. ||i,j and in Table 0). 

Let us now consider the metallicity gradients in the 
subsystems. Fig. || displays the R a — [Fe/H] and Z max — 
[Fe/H] diagrams. The straight lines are the least-square 
regressions. Assuming that the stars are born near the 
apogalactic radii of their orbits, the slopes of these lines re¬ 
flect the initial radial and vertical metallicity gradients for 
the subsystems. To increase the accuracy of the gradient 
estimates, we rejected the most distant data points in each 
case (see the vertical dotted lines in the diagrams). The 
resulting gradients are presented in Table 1 and the corre¬ 
lation coefficients in the corresponding diagrams. Only the 
vertical gradient in the proto-disk halo exceeds the error, 
but the correlation coefficient given in Fig. |]c tends to zero 
if we add the rejected distant stars. Note that RR Lyrae 
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Fig. 2. Distributions of orbital elements for the RR Lyrae stars. Result for stars with V^, e c < 280 kms 1 (left 
column), V pec > 280 (right colimn). The curves in the top row represent an approximation of the distributions of 
Gaussian functions. In bottom row, the dots are the reconstructed distributions of RR Lyrae stars in Z (see the text 
for details), and the curves approximate the reconstructed distributions with an exponential law (the corresponding 
scale heights and their uncertainties are indicated). 
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Table 1. Characteristics of metal-poor subsystems of field RR Lyrae variables and Globular Clusters 


Characteristics 

Proto-disk halo 

Accreted halo | 


RR Lyr 

GC 

RR Lyr 

GC 

< [Fe/H] >, dex 

-1.56± 0.04 

-1.71i 0.05 

-1.62i 0.03 

•4:60= 0.07 

V[Fe/H], dex 

0.33± 0.05 

0.26i 0.03 

0.32i 0.02 

0.35i 0.05 

< 0 >, km/s 

93± 8 

77± 33 

-51i 11 

-23i 54 

ere, km/s 

83± 6 

129± 19 

106i 8 

140i 18 

< e > 

0.57±0.03 

0.53i0.06 

0.76i0.02 

0.59i0.06 

limRa, kpc 

18 

10 

42 

20 

lim Z max , kpc 

9 

10 

22 

20 

Zo, kpc 

2.6i0.2 * 

2.5i0.5 

6.1i0.3 

8.5il.5 

gradit[Fe/H], kpc^ 1 

-0.00i0.01 

-0.03i0.02 

-0.00i0.01 

-0.03i0.01 

grad z [Fe/H], kpc -1 

-0.02i0.01 

-0.03i0.02 

-0.00i0.01 

-0.03i0.01 


Note.* - The gradient of the proto-disk halo is strongly distorted by the selection which depends on only nearest stars. 


Proto-disk halo 


Accreted halo 



Z 


max 


, kpc 


Fig. 3. Metallicities of the RR Lyrae stars versus apogalactic distance (a,b,) and maximum height above the galactic 
plane (c,d). The straight lines are least-square fits. The slopes of the lines determine the metallicity gradients. The 
most distant stars in each panel (on the right of the dotted lines) were rejected from our computations. The correlation 
coefficient and their uncertainties are given. 

metallicities from Zinn & West (1984) are estimated using al. 1997, Rey et al. 2000). Therefore, to check our result 
the AS method, which actually measures calcium abun- we have to use more reliable metallicity determinations for 
dance. According to current evidence there is a system- field RR Lyrae stars, when sufficient data will be available, 
atic uncertainty between the metallicity scale suggested Thus the existence of vertical gradient in the proto-disk 
by Zinn & West (1984) and other ones (see Rutledge et halo is still an open question. In contrast, according to our 
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investigation, the complete absence of a radial gradient in 
the proto-disk halo and of both gradients in the accreted 
halo for the RR Lyrae stars is beyond doubt. 

Suntzeff et al. (1991) investigate radial and vertical 
nretallicity gradients for the metal-poor field RR Lyrae 
population as a whole. They used the present positions 
of stars in selected regions of the Galaxy in the galac- 
tocentric distance range 4-30 kpc. The authors indicate 
that outside the solar circle the metallicity gradients are 
zero as a function of R a or \Z\. This result is in agree¬ 
ment with gradients in our accreted halo estimated on R a 
and Z max . Inside the solar circle they find the radial gra¬ 
dient of —0.06 ± 0.02 dex kpc -1 , which is steeper than 
that of our proto-disk halo for field RR Lyraes (but is in 
agreement for globular clusters). This discrepancy is due 
to the lack of observed RR Lyraes inside the solar circle. 
The vertical gradient in [Fe/H] they estimate as roughly 
-0.05 dex kpc -1 out to 3 kpc. Our value for the proto¬ 
disk halo is larger because Suntzeff et al. (1991) did not 
divide the metal-poor stars into subsystems. So we may 
conclude that computed orbital elements for the nearest 
stars are good enough for estimation of metallicity gradi¬ 
ents in the much larger distance range (but not inside the 
solar circle.) 

5. Discussion and conclusion 

Here, we will compare the characteristics of the metal-poor 
RR Lyrae subsystems derived in this study to the parame¬ 
ters of corresponding subsystems of globular clusters from 
Borkova & Marsakov (2000), since only globular clusters 
were distinguished according to an intrinsic, physical pa¬ 
rameter rather than interrelated spatial and kinematical 
criteria. The parameters of the metallicity distributions 
of the corresponding halo subsystems differ somewhat. In 
particular, the mean metallicity of the proto-disk halo de¬ 
rived from the globular clusters is lower than that of the 
outer halo. The metallicity dispersion is also lower. The 
field stars show the opposite pattern (Table |l]) . In all cases, 
however, the differences are comparable to the formally 
computed uncertainties, indicating that any conclusions 
about differences between these parameters have low sta¬ 
tistical significance. The vertical gradients in the proto¬ 
disk halo for both the RR Lyrae stars and the globular 
clusters almost coincide, but the radial gradients differ 
(see Table [[]). In the accreted halo, both gradients are ab¬ 
sent for the field RR Lyrae stars but are non-zero for the 
clusters. However, both gradients for the globular clusters 
in the accreted halo are due exclusively to metal-richer 
objects close to the galactic center (R ~ 7 kpc). Distant 
RR Lyrae were not included in our sample. In any case, 
values of all corresponding gradients coincide within the 
uncertainties. In this study, we have identified the halo 
subsystems based on spatial velocity. Therefore, the differ¬ 
ences between the proto-disk halo and the accreted halo in 
any kinematical parameter for field stars should be more 
prominent. Indeed, while the difference between the or¬ 
bital velocities for the globular clusters of the proto-disk 


halo and of the accreted halo is ~ 100 kms -1 this differ¬ 
ence is approximately 40% higher for the RR Lyrae vari¬ 
ables (see Table||). The velocity dispersions for the sub¬ 
systems of globular clusters are obviously overestimated 
due to the large distance uncertainties and, as a result, 
are much higher than the values for field RR Lyrae stars. 
The mean eccentricities in the proto-disk halo subsystems 
are the same, whereas in the outer halo the eccentricities 
are, on average, higher for the field stars, as expected (note 
that the proper motions and, hence, orbital eccentricities 
are known only for a small number of clusters, and with 
large uncertainties). The radial size of the proto-disk halo 
is approximately a factor of 1.8 larger for the field stars 
than for the clusters, whereas the two scale heights were 
the same within the errors. Recall that we can estimate 
the radial sizes of RR Lyrae subsystems only from their 
maximum distances from the galactic center, which leads 
to appreciable overestimates of these sizes. The radial and 
vertical sizes of the outer accreted halo subsystem of field 
stars are naturally the largest, and are in reasonable agree¬ 
ment with the corresponding sizes for the subsystem of 
globular clusters. Note that, in order to obtain correct es¬ 
timates of sizes of galactic subsystems based on data for 
nearby stars, it is necessary to take into account the kine¬ 
matical selection effect, which leads to a deficiency of stars 
with large R a and Z ma x in the solar neighborhood. 

Thus, we found that the generally good agreement be¬ 
tween the characteristics of corresponding subsystems of 
field RR Lyrae stars and globular clusters, distinguished 
using different criteria, shows that both populations are 
not uniform. Both the clusters and field stars belong to 
two spherical subsystems of the Galaxy: the inner proto¬ 
disk halo (related to the disk by its origin) and the outer 
accreted halo. The collected results indicate that this sub¬ 
system is characterized by large size, an absence of appre¬ 
ciable metallicity gradients, predominantly large orbital 
eccentricities, a large number of objects on retrograde or¬ 
bits, and, on average, younger ages for its objects, support¬ 
ing the hypothesis that objects in this subsystem have an 
extragalactic origin. 
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